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Abstract. The kinematics and physical conditions of 
the bipolar planetary nebula M 2-48 are analysed from 
high and low dispersion long-slit spectra. Previous CCD 
narrow-band optical observations have suggested that this 
nebula is mainly formed by a pair of symmetric bow- 
shocks, an off-center semi-circular shell, and an inter- 
nal bipolar structure. The bipolar outflow has a complex 
structure, characterised by a series of shocked regions lo- 
cated between the bright core and the polar tips. There 
is an apparent kinematic discontinuity between the bright 
C^Q [ bipolar core and the outer regions. The fragmented ring 

■ around the bright bipolar region presents a low expansion 
[ velocity and could be associated to ejection in the AGB- 

■ PN transition phase, although its nature remains unclear. 
O ' The chemical abundances of the central region are derived, 

■ showing that M 2-48 is a Type I planetary nebula (PN) . 



keywords: ISM: jets and outflows — echelle spec- 
troscopy — Hydrodynamics — Shock waves — planetary 
nebulae : individual: M 2-48 
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ASTRONOMY 

AND 
ASTROPHYSICS 



1. Introduction 

The morphology of M 2-48 has been previously dis- 
cussed by Vazquez et al. (2000) [hereafter Paper I] who 
found the full extent of its bipolar structure and indica- 
tions of the presence of high-velocity outflows from their 
monochromatic images. The occurrence of highly colli- 
matcd, high velocity jet-like components in planetary neb- 
ulae is well established (e.g. Miranda & Solf 1992; Lopez 
et al. 1995; Sahai & Trauger 1998; Miranda, Guerrero & 
Torrelles 1999; O'Connor et al. 2000; Miranda, Guerrero & 
Torrelles 2001 and Gongalves, Corradi & Mampaso 2001, 
etc.). These bipolar outflows usually show indications of 
episodic outburst and rotation or precession of the symme- 
try axis (e.g. Lopez et al. 1997). The origin of such outflows 
is intimately related to the formation and shaping of PNe 
and thus their kinematic characterisation is a key param- 
eter for the models. Several models have been proposed 
for the generation of collimated outflows in PNe, such as 
pure hydrodynamical models (Icke et al. 1992, Mellema 
1995, 1997), models involving stellar rotation and mag- 
netized winds (Rozyczka & Franco 1996; Garcfa-Segura 
1997; Garcia-Segura & Lopez 2000), and the influence of 
binary system (Soker & Livio 1994; Reyes-Ruiz & Lopez 
1999, Soker 2002). 
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The bipolar planetary nebula M 2-48 has three main 
structural components: (1) an internal bipolar central re- 
gion (CR); (2) an off-center semi-circular fragmented shell, 
containing a region with a bow-shock morphology (this 
has been interpreted in Paper I as the possible interac- 
tion of the bipolar, collimated outflow with the shell), and 
(3) a pair of symmetric knots with a seemingly bow-shock 
morphology, located at either side at « 2' from the center 
(see Paper I). 

In order to analyze the kinematical structure, exci- 
tation mechanisms and chemical abundances of M 2-48, 
spatially resolved high and low dispersion spectroscopy 
have been now obtained. We first describe the observa- 
tions and the results extracted from the position-velocity 
(PV) diagrams in section |^. The kinematical structure and 
morphology are discussed in section ^, including a simple 
model to predict the orientation of the shocked structures 
with respect to the plane of the sky. In section fl|, the ex- 
citation mechanism and chemical abundances at different 
regions of the nebula are derived. Finally, our main con- 
clusions are summarized in section ^. 

2. Observations and Results 

Low-dispersion (LD) spectroscopy was obtained in 1999, 
July 18, 19 and 20, with the B&Ch spectrometer com- 
bined with the f/7.5 focus of the 2.1-m OAN telescope 
at the Sierra de San Pedro Martir, B. C, Mexico. A 300 
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Fig. 1. Slit positions of the low-dispersion spectroscopy Fig. 3. Slit positions of the high-dispersion spectroscopy 
are marked over the contour plot of a [N II] 6584 A CCD (A-F) are marked against a contour plot of a [N II] 6584 A 
image of M 2-48 taken from Vazquez et al. (2000). CCD image of M 2-48 taken from Vazquez et al. (2000). 



Table 1. Data of the ID spectra extracted from the low 
dispersion spectra. 




5000 6000 

Wavelength (A) 



Fig. 2. Low-dispersion optical spectra of two different re- 
gions of M 2-48. The central region (CR) is shown ex- 
panded along the flux axis in the middle panel. The bow- 
shock Bl is shown in the bottom panel. 



lines/mm grating was combined with a 220 /im (= 2.9 arc- 
sec) slit width and a Tektronix CCD of 1024 x 1024 pixel, 
24 /im square pixels (= 0.84 arcsec) as detector. The re- 
sulting spectra cover from 3400 to 7500 A, with a 12 A 
spectral resolution. The slit was oriented in several posi- 
tion angles (PA) to cover different regions of the PN. Fig- 
ure [l] shows the slit positions over a contour plot of the 
[N II] 6584 A image (Paper I). The LD CCD frames were 
trimmed, bias-subtracted, flat-fielded and sky-subtracted 
with standard techniques using IRAF. Flux calibration 
was performed using sensitivity functions derived from the 
standard star HD 192281. Table |l| presents the log for each 
of the ID spectra that have been extracted and analyzed 
and Figure || shows the one-dimensional spectra of CR 
and Bl. 





RA. 


Exposure 


Slit length 


Region 


n 


time (s) 


(arcsec) 


CR 





3600 


25.9 


si 


32 


2400 


11.5 


s4 


-32 


1800 


8.6 


Bl 


47 


1800 


20.2 


Bl 


50 


3000 


4.8 



The high-dispersion (HD) spectroscopy was obtained 
in 1999, June 29 and 30, with the Manchester Echelle 
Spectrometer (MES; Meaburn et al. 1984) combined with 
the f/7.9 focus at the 2.1-m OAN telescope at the Sierra 
de San Pedro Martir, B. C, Mexico. This spectrometer 
has no cross-dispersion. A 90 A bandwith filter was used 
in order to isolate the 87th order, which contains the Ha 
and [N II] 6548, 6584 A emission lines. 

A Tektronix CCD with 1024 x 1024, 24 [im (= 0.3 arc- 
sec) square pixels was the detector. A 2 x 2 binning was 
employed in both spatial and spectral dimensions. The 
"seeing" varied between 1 and 2 arcsec during the ob- 
servations. In Figure 0, the slit positions A-F are shown 
against a contour plot of [N II] 6584 A of M 2-48. The 
150 fim (= 1.9 arcsec) wide slit was oriented at a position 
angle PA=75° for slits A - D and North-South for slits E 
-F. The spectral resolution is lOkms -1 . Exposure times 
of 1800 s were carried out for each slit position. Images 
were reduced in the standard manner using IRAF. The 
spectra were wavelength calibrated with a Th-Ar lamp to 
an accuracy of ±1 kms -1 . 

The bidimensional arrays of [N II] 6584 A position- 
velocity (PV) line profiles from the HD spectra (slits A-D), 
are shown in Figure |^. Figures 5, 6 and 7 show enlarge- 
ments of the core region for slits B and C; the line profiles 
from slits E and F and the line profiles at the tip of north 
eastern lobe for slits B and C, respectively. 
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Fig. 4. Gray-scale representation of the position-velocity 
arrays of [N II] 6584 A along slits A - D. Labels correspond 
to the regions indicated in Figure 3. 



3. Morphology and Kinematics 

M 2-48 has a bipolar morphology, as has been described 
in Paper I (see also the images in Manchado et al. 1996). 
From a bright bipolar core emerge relatively faint, bro- 
ken filaments that delineate fairly tight collimated lobes 
extending over 2 arcmin on either side. Some of the fila- 
ments close to the core seem to trace an incomplete arc 
that gives the appearence of being the remnants of a spher- 
ical shell. The PV diagrams shown in Figures 4-7 reveal 
the complex kinematic structure of the main components 
of M 2-48. The labels Bl, si, s2, s4 and s5 in the PV ar- 
rays in Figure |J match the location of the regions with 
the same labels in Figure 3, where we have followed the 
nomenclature of Paper I. The spectral element labeled as 
(fl) is probably related to the tip of the southwestern lobe. 

3.1. Central region 

Line profiles from the bright nebular core are shown in 
detail in Figure ^. The systemic velocity of the core, as 
defined by the central intensity maximum of the line pro- 
files shown in Figure |^ lies around —15 km s _1 . The line 
profiles replicate a bipolar outflow emerging from a bright, 
compact core, extending bluewards to -50 km s _1 on the 




Fig. 5. [N II] 6584 A PV dia grams of the central region 
of M 2-48 (for slit positions B and C). The line profiles 
are shown with a greyscale and overlayed logarithmic con- 
tours. 




Fig. 6. [N II] 6584 A PV dia grams of the bipolar lobes of 
M 2-48 (for slit positions E and F). The line profiles are 
shown with a greyscale and with logarithmic contours. 



west side and redwards to +25 km s _1 on the east side. 
The line profile from slit C (bottom panel in Fig. 5) shows 
particularly extended wings, reaching remarkable values 
of ±85 km s -1 at FWZI within the central (~ 20 arcsec) 
region. The bright compact core itself is seen elongated 
and tilted in the velocity space, expanding towards nega- 
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Fig. 7. [N II] 6584 A PV diagrams of the bow-shock Bl 
of M 2-48 (for slit positions B and C). The line profiles are 
shown with a greyscale and with logarithmic contours. 



tive velocities on the east side (—25 km s 1 ) and positive 
velocities on the west side (+10 km s _1 ). 

3.2. The arc-like structure 

One of the curious features of M 2-48 is a group of 
several knots around the central bipolar region that seem 
to form an arc-like structure (Paper I). These knots (sl-s5 
in Figure |IJ) were interpreted as possible segments of an 
expanding semi-circular shell. Some slits were obtained 
with the purpose of intersecting this possible expanding 
shell and clarify their nature. The line profiles shown in 
Figure are found to follow the velocity pattern of the 
bipolar outflow emerging from the extended regions of the 
core, with redshifted velocities on the east side (« +10 km 
s _1 ) and blueshifted velocities on the west side (« —25 km 
s _1 ). No indications of a velocity ellipse that could hint to 
the presence of an expanding shell are detected, at least 
with the present spectral resolution. 

3.3. Bow-shock si 

The eastern lobe is the brighter one, with several high- 
speed knots with bow-shock looking morphologies being 
detected in the HD spectra, namely Bl, si and s2. On the 
western lobe we detect and identify knots s4 and s5, and 
the feature f 1 (see Figures ||| and ^) . 



A high-velocity element, labeled si in Figure 4, coin- 
cident with the position of a bright knot on the eastern 
side and near a portion of the arc-like structure, shows an 
apparent bow- shock morphology and kincmatical struc- 
ture (see Figures [| and ^ bottom right and the image in 
Manchado et al. 1996). It is of interest to estimate the 
orientation angle of the bow-shock associated with knot 
si with respect to the plane of the sky, <f>. To do that, 
we have used a "3/2-D" model to predict the position- 
velocity diagram of a bow-shock at different angles with 
respect to the plane of the sky. Raga & Noriega- Crespo 
(1993) found that the Ha emission per unit time and area 
is approximately given by : 

Iu a oc (v shock ) 7Ha , (1) 
where v s hock is the shock velocity and 7H a is a constant 
dependent on the shock velocity regime. 

If we assume that the jet is moving along the z-axis, 
the shape of the bow-shock can be parametrized by : 
r«z p , (2) 
where p determinates the aperture of the bow-shock. We 
have used a value of p — 2. 

Thus, the input parameter for the model essentially is 
the shock velocity. This velocity can be estimated from 
the PV diagrams obtained from the long-slit spectroscopy 
(Figure |). 

In addition, we use a value of 7h q = 3.49 correspond- 
ing to the 20.0-81.3 km s~ x range of velocities proposed 
by Raga & Noriega-Crespo (1993). 

We can calculate the Ha emission, and predict the 
position- velocity diagram at different orientations with re- 
spect to the sky plane, as shown in Figure |[ From this 
comparison we can establish that si is inclined at an angle 
^10° with respect to the plane of the sky. 

This low inclination angle indicates that the radial ve- 
locity of the outflow producing the si bow-shock can be 
a factor of about 6 larger than the velocities derived from 
the PV diagrams. 

We have used the radial velocity of the peak emission 
of si, 15 km s _1 , as the expansion velocity projected 
to the line of sight. Using the inclination estimated for 
si, <fi ~ 10°, we can determine the tangential expansion 
velocity. Finally, taking into account the projected linear 
separation between si and the central region we can esti- 
mate a dynamical time of r s i « 3 [d/kpc] x 10 3 yr. 

In addition, a peculiar feature labelled f2 can be dis- 
tinguished in the P V diagram of slit C (Figure [|) , with a 
constant heliocentric velocity of about —60 km s -1 . This 
feature seems to be spatially coincident with the eastern 
arc, corresponding to the diffuse material between knots 
si and s2. 

3.4- External bow- shock Bl 

In the [N II] 6584 A direct images (see Paper I), it is 
possible to see two symmetrical structures with a bow- 
shock morphology. The eastern knot, Bl, which is brighter 
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Fig. 8. Comparison of the position- velocity diagram of the 
si bow-shock and "3/2-D" bowshock model fits at differ- 
ent orientations pointing away from the plane of the sky. 
The best fit is at = 10° . 
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Fig. 9. Comparison of the position- velocity diagram of the 
Bl bow-shock and "3/2-D" bowshock model fits at differ- 
ent orientations pointing inwards to the plane of the sky. 
The best fit is at = 10°. 



than the western one, B2, is blueshifted and has a line 
width of w 80 km s^ 1 . The PV dia grams of this bow-shock 
are shown in Figure where we can see two different zones 
of Bl covered by the slits B and C. Slit B slightly covers 
the center of the bow-shock (larger line width), and C 
covers the wings (smaller line width). 
It is possible to use the "3/2-D" model described in section 
3.2 to predict the PV diagram of this bow-shock. From the 
comparison of the observations and the numerical model 
of the bow-shock Bl for different orientations (see Fig- 
ure ||), we can estimate that Bl is inclined at an angle of 
~ 10° pointing inwards to the plane of the sky. The angle 
with respect to the plane of the sky suggests that the ra- 
dial velocity of the bow-shock Bl can be a factor of w 6 
greater than the velocities derived from the PV diagrams. 
In order to establish a sequence of e vents, using the same 
procedure as for si (see section |3~3|) , we estimate the dy- 
namical age of this bow-shock to be tbi ~ 5.3 [d/ kpc] x 10 3 
yr. This is consistent with the idea that the bow-shocks 
Bl and si correspond to different bipolar ejections and 
probably with different axis orientation. This suggests a 
change in the orientation of the symmetry axis of the suc- 
cesive outflows. In this way, a previous ejection of material 
would form the expanding arcs interacting with a bipolar 
outflow and the most external bow-shocks Bl and B2, and 
a later high-velocity bipolar outflow with different orien- 
tation would form the si bow-shock. 



We have also noted that the symmetry axis of the 
most internal bipolar region is not coincident with either 
of these two high- velocity bipolar outflows. 

In addition to the main features described in previous 
sections, we can see other interesting features in the PV 
diagrams obtained with the HD spectroscopy. In Figure |^ 
we can distinguish a very weak feature labelled fl in slits B 
and C. This is a component with an heliocentric velocity of 
« —50 km s . This emission could correspond to the wing 
of the western bipolar lobe outflow. Observations at other 
slit positions or longer exposure times would be necessary 
to resolve in more detail the kinematical behavior of this 
feature. 



4. Line intensity ratios and chemical abundances 

As noted in section ||, five regions were selected 
from the slit positions observed with low resolution spec- 
troscopy (see Figure Q), covering the four different mor- 
phological zones pointed out in Paper I. In Figure |^, the 
flux calibrated spectra of two different zones are shown. 
The top and middle panels correspond to different scales 
(high and low intensities, respectively) of the spectrum of 
the central region. The bottom panel corresponds to the 
spectrum of Bl (PA 47°). 

The extinction coefficient, C(H/3), was derived for 
three regions: Central Region (CR), si and Bl (PA 47°) 
from the comparison of observed and theoretical Ha/H/3 
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Table 2. Dereddened line intensity ratios. 
A - Bright regions 
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ratios. We have used the calculations of Storey & Hum- 
mer (1995) assuming case B and an iteration procedure 
adopting finally the electron densities (iV e ) derived for 
each particular zone and the electron temperature (T e ) 
obtained for CR (T e =10800 K) for all of the regions. The 
reddening law of Seaton (1979) has been used. Dereddened 




0.5 1 : 

log [Ha /[SII] (6717 + 6731)] 



Fig. 10. Diagnostic diagrams showing the emission line 
ratios of the different regions observed in M 2-48. Top: log 
Ha/ [SII] vs. log Ha/ [Nil]. Bottom: log Ha/ [SII] vs. the 
electron density indicator [SII] 6717/6731 (adapted from 
Sabbadin et al. 1977) 



Table 3. Chemical abundances in the central region 



12+log 0+/H+ 


7.85 


12+log 0++/H+ 


8.43 


12+log O/H 


8.55 


log N+/0+ 


+0.30 


12+log Ne++/H+ 


7.75 


12+log S+/+ 


6.16 


12+log Ar ++ /H+ 


6.30 


12+log Ar 3+ /H+ 


5.77 


He+/H+ (4471) 


0.127 


He+/H+ (5876) 


0.098 


He+/H+ (6678) 


0.096 


He+/H+ (7065) 


0.163 


<He+/H+> 


0.100 


He ++ /H+ 


0.008 


He/H 


0.108 



line intensity ratios with respect to H/3 as well as the val- 
ues of C(H/3) for these regions are presented in Table |[ 
The uncertainties in the line fluxes are of the order of 5% 
for the brighter lines and about 20-30% for the weaker 
ones. Emission line fluxes of blended lines such as the 
Ha+[N II] nebular lines and the [S II] doublet have been 
obtained from Gaussian fitting of the line profiles using 
the Starlink DIPSO package (Howarth & Murray 1990). 
The dereddened, integrated H/3 flux of each region is also 
included in Table 0. 
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log j[SII] (6717 + 6731)/Haj 



Fig. 11. Diagnostic diagram representing log [SII]/Ha vs. 
log [OIII]/Ha of the different regions observed in M 2-48 
(adapted from Phillips & Cuesta 1999) 

N e and T e (see Table ||) have been derived following 
standard methods of collisionally excited line ratios us- 
ing the five-level program for the analysis of emission-line 
nebulae of Shaw & Dufour (1995). There seems to be a 
gradient in N e , with parameter decreasing values towards 
the external zones. This is consistent with the expansion 
of stellar ejecta as the distance towards the nucleus in- 
creases. T e derived from [O III] and [N II] lines for CR 
are almost coincident, indicating that there is no strong 
temperature stratification in the bright bipolar nebular 
core of the nebula. The T e derived for si seems too high 
to be accounted for by pure photoionization, confirming 
contamination by shock emission in this zone, which mor- 
phologically and kincmatically resembles a bow-shock (see 
section 



3.3) 



The spectra of the other two zones: s4 and Bl (PA 
50°), are so faint that H/3 is not measured and therefore it 
is not possible to derive the extinction coefficient for these 
regions. The emission line ratios with respect H/3 and the 
uncorrected Ha flux for these regions are included in Table 

We have used several diagnostic diagrams in order 
to evaluate which is the dominant excitation mechanism 
prevalent in these regions. In Figure [l^ we show two dia- 
grams involving the locus of several line ratios: Ha/[N II], 
Ha/[S II], and [S II] 6717/6731 for H II regions, PNe and 
supernova remnants (adapted from Sabbadin et al. 1977). 
These diagrams clearly indicate that CR is mainly radia- 
tively excited. However, the positions of si and Bl are not 
so definite, in the sense that some shock contribution is 



present. The fact that the positions of Bl and si are just 
below the locus of the SNRs in the upper pannel of Fig- 
ure [l(] could be also due to the high N enrichment of the 
nebula (see below). 

In Figure |ll| we show a diagram adapted from Phillips 
& Cuesta (1999) involving [O III] /Ha and [S II] /Ha. This 
diagram includes the locus of the observed emission line 
ratios for a large number of radiatively excited PNe and 
predictions from plane-parallel and bow-shock models by 
Hartigan et al. (1987) and Shull & McKee (1979). The 
position of the different regions of M 2-48 in this diagram 
confirms that CR is radiatively excited and that shock 
excitation is contributing to the spectra of si and Bl. As- 
suming purely shock excitation for these zones and com- 
paring their line intensities with the models by Hartigan 
et al. (1987) we find the best agreement for models with 
shock velocities of 60 and 130 km s _1 , respectively, in 
agreement with the values obtained in previous sections. 

Ionic abundances of the CR are shown in Table |3| and 
have been derived using the emission line ratios given in 
Table ||, the code of Shaw & Dufour (1995), and assum- 
ing T e = 10 800 K and N e = 1260 cm" 3 . To derive the 
He + /H + and He ++ /H + ratios we have used the effective 
recombination coefficients given by Pequignot et al. (1991) 
for all the lines involved. The He + abundances obtained 
from each individual line have been corrected for colli- 
sional contributions following Kingdon & Ferland (1995). 
The final, adopted value of the Hc + abundance is the av- 
erage of the values obtained for the different individual 
lines excluding Hel 7065 A, which suffers from the largest 
collisional effects. 

The total O abundance has been determined from the 
derived abundance of + and ++ combined with the 
ionization correction factor (ICF) for the presence of 3+ 
proposed by Kingsburgh & Barlow (1994). The total He/H 
ratio has been obtained simply by adding the <He + /H+ > 
and Hc ++ /H + ratios. 

The log N+/O+=+0.30 obtained for the nebula clearly 
indicates that it is a nitrogen-enriched object. Taking into 
account the chemical classification proposed by Peimbert 
& Torres- Peimbert (1987), M 2-48 should be classified as 
a Type I PN (log N/O > -0.30 and/or He/H > 0.125), 
although it is not a helium-rich object. The O/H ratio is 
slightly lower than the average of the Galactic Type I PNe 
(12+log O/H = 8.65, Kingsburgh & Barlow 1994). 

5. Conclusions 

The kinematic mapping of M 2-48 has revealed a com- 
plex bipolar structure characterised by a series of shocked 
regions extending accross the major axis of the nebula 
and reaching the tips found in paper I. Indications that 
the bipolar outflow has suffered some changes in direction 
have also been detected. 3/2 D shock models of the line 
profiles are consistent with the main bipolar outflow hav- 
ing a small ± 10° angle with respect to the plane of the 
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sky, thus making the real outflow velocity some six times 
larger than the observed one. The line fluxes derived from 
the low dispersion spectrum in the central region show 
that M 2-48 is nitrogen enriched; and should be classified 
as type I. 
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